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Abstract 

We propose and develop a model as an attempt to describe the 
local mechanism of mechanical alloying. This model is based on the 
observation of global parameters, such as the volumes of the various 
phases present in the material, and their transformations during the 
process. The model is applied to milled Fe-Cu and Fe-N experimental 
results obtained in previous works ^|2j. In the milling of Fe-Cu sys- 
tem, four stages, whose features are in accordance with experimental 
observation, are identified: (i) energy storage; (ii) initiation of reaction 
and rate increase; (iii) decrease of the reaction rate and (iv) station- 
ary stage. No energy storage stage is observed and no saturation 
is achieved in the milling of Fe-N system. Moreover, we propose a 
modification on the global kinetic law introduced in [T2| based on the 
different nature of Fe-Cu and Fe-N systems. This modification allows 
the extension of the validity of this law. 
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1 Introduction 



In classical studies of phase transformations, kinetic laws allow a global de- 
scription of volume phase as a function of time evolution thus distinguishing 
them from growth laws. However, growth laws translate a local microscopic 
modeling of phase transformation phenomena and provide the construction of 
kinetic laws from microscopic observation and (or) models. Both approaches 
exist in mechanical alloying (MA) although they are rarely explicit in the 
literature. 

MA is a cheap, efficient, and flexible technique (see the survey by Surya- 
narayana j3]) which has attracted much attention from scientists of various 
research areas. Recently, MA has been used to produce materials for biolog- 
ical application j4j and to investigate the spin-glass character of nanograin 
boundaries [H], amongst many other applications. Although much of the re- 
search is based on empirical conclusions obtained directly from experimental 
results, there has been some attempts to make a quantitative description of 
MA. A quantitative description is very important and necessary in order to 
provide a way to establish predictable capabilities for the MA processes 

Abdellaoui and Gaffet [71IH] have shown that neither the shock energy nor 
the shock frequency govern the NiioZry amorphization process in a planetary 
ball mill if taken into account separately. Only the injected shock power 
is responsible for the ball milled final product. In another study of Ni-Zr 
amorphization, Chen et al. using a vertical vibrating frame grinder, have 
found that these compounds become fully amorphized if a defined milling 
intensity (/) is greater than a threshold value. The milling intensity is defined 
as: 

" M, ' 

where is the ball mass, Vmax is the maximum velocity of the ball, / the 
impact frequency, and Mp the mass of the powder in the vial. Therefore, 
according to [9' , the same final product can be achieved by transfering large 
momentum with low frequencies or small momentum with high frequencies. 
These two approaches, while successful in describing individual processes in 
particular types of mill, have some built in difficulties that prevent them to 
fully describe a general MA experiment. The shock power proposed in [71IH] 
is a function of balls' kinetic energy and impact frequency. However, in a 
planetary mill, both parameters are functions of the velocity of the balls in 
the sense that the faster a ball moves, the higher its kinetic energy and the 
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more times it will suffer impacts. Moreover, in the vertical vibrating mill 
used to investigate the role of the milling intensity [HI, the shock frequency 
is determined by the inverse of the time taken by ball to move up and down. 
This time is proportional (and as a consequence the frequency is inversely 
proportional) to potential energy. The velocity Vmax is itself proportional to 
potential energy which makes the product Vmaxf close to a constant value. 
These two approaches also neglect the importance of the material's response 
to the impacts. 

MA is a "far from equilibrium technique" and despite its dynamic char- 
acteristics Yavari et al. ^0, JJJ have employed equilibrium thermodynamics 
to explain the MA mechanism. Pure elements form nanocrystalline struc- 
tures after severe mechanical deformation and the energy stored in the large 
grain boundaries could be the driving force during alloying. Examining the 
immiscible Fe-Cu system prepared by mechanical alloying, Yavari et al. have 
shown that the available interfacial energy is not enough to compensate the 
positive enthalpy of mixing. However, repeated fractures break the particles 
into very small fragments with high free energies which dissolve themselves 
into the matrix by the Gibbs-Thomson effect. 

Matteazzi et al. ^2 have proposed a global kinetic law for MA of iron 
carbides: 

V{t) = ke-''^ (2) 

where V{t) is the transformed phase volume and k, b and c are empirical 
constants. They have proposed that c is connected to impact energies and b 
is related to "work conditions" . The value of k sets the scale. 

In a previous work P, we reported a discussion about the MA pro- 
cess based on experimental and numerical results. A nano cry st aline fcc- 
FeCu solid solution was obtained by MA and analyzed by x-ray diffraction, 
Mossbauer spectroscopy, and differential calorimetry. In addition, simula- 
tions of the material behavior during uniaxial deformations, compatible with 
what we expect on an impact were performed using molecular dynamics. We 
claimed that the number of balls in the mill vial, which is related to impact 
frequency, did not affect the energy transfer processes, but rather only the 
reaction rates. This finding has recently been confirmed by Girot et al. [TB] . 
We also showed that a local description of MA process is necessary in or- 
der to understand its driving mechanisms and that many factors such as; 
powder deformation by the milling tools, chemical and physical properties 
of materials, strain/stress fields, and collision energies have to be taken in 



3 



account together to describe MA. This view was recently raised in the review 
by Takacs [T^ . 

In this contribution we propose and develop a model as an attempt to 
describe the local mechanism of MA. This model is based on the phase volume 
transformations during the process and is applied to milled Fe-Cu and Fe-N 
experimental results obtained previously [I1I21- 



2 The Model 

Our first step to understand MA mechanisms is to investigate how the phase 
transformations are driven in the course of the process. This can be done 
using measured global parameters such as the phase volumes. If a given phase 
a occupies a volume V^(t) in the total volume V ^ then the concentration 
Xait) is given by: 

^.W = (3) 

MA is essentially a discrete process characterized by discrete energy trans- 
fers taking place at individual impacts ^21, which means that only a small 
fraction of the powder is processed at each event. Moreover, a particle of 
powder under an impact has no memory from the previous impacts Sig- 
nificant powder transformations occurring during milling impacts are con- 
fined in a small volume v which is the powder volume trapped by the milling 
tools, as suggested in [12]. The a phase volume in v is and assuming a 
homogeneous powder JH], the relative concentration is then given by: 

X^{t) = (4) 

Let us consider that a phase is created after a MA driven solid state reac- 
tion between two other phases f3 and 7. The volume created or transformed 
Ava after an impact involving a jSj interface is a function of the probability 
to find a phases reaction (Pp^it)) and tribochemical activity (^^^(t)) given 
by: 

Av.it) ^F{P^,it),Cp^{t)). (5) 

The word tribochemical includes many aspects such as chemical activities, 
morphology of components and impact characteristics (energy, frequency. 
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etc.). In a first approximation eq. |31can be written as: 

Ava{t) = J2Pp,mUt)v. (6) 

13-y 

This summation is carried out over all f3'j phase reactions which could pos- 
sibly lead to production of a. 

A formal description of ^p^(t) is quite difficult; however, it is possible 
to describe the volume of a phase during the total milling time At. The 
evolution of the a phase global volume (AV^(t)) in the course of time interval 
At is given by: 

AV^it) = Av^it) X fAt = Y^Pp,{t)Cp^{t)vfAt, (7) 

Pi 

where / is the impact frequency and f At is the number of impacts during 
the time interval At. 

The probability to find a /?7 interface which can react is given by the 
relative surface areas of the j3 and 7 {Sj3{t) and Sy{t)) phase grains in the 
volume v: 

P,.(t) = ^i^^^, (8) 

where 5*^ is the total surface area in the volume v. Considering spherical 
grains, eq. |H1 can be written as: 

where rj is the particle radius of phase i and the summation is carried out 
over all the phases in the volume v. 

This analysis can be greatly simplified if we consider two-phased systems, 
composed by phases e and 6. The possible reactions between the phases are 

e + 6 — y Ave = 0, Avs = 0, (10) 

S + 6 — y Ave = 0,Av5 = 0, (11) 

e + 6 — > Av^ <0, Avs > or Av^ > 0, Avs < 0. (12) 

The summation in eq. [Z|is reduced to just one term to translate the e (or 6) 
volume phase evolution in such a way that: 

AK(t) = Pesmisit^f^t- (13) 
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In this way, the tribochemical activity Ctsit) can be written as: 

_ AK(t) 1 

At P,s{t)vf ^^^^ 

Equation IT^ presents at the same time a kinetic factor ( ^ a/ ) ' ^ term of 
chemical nature {P^sit)) and a "condition of milhng" term (vf). The prob- 
abihty of a e6 interface is given by eq. IHl and can be written as: 

and if the particles radii are comparable (r^ ^ r^), eq. HSlcan be reduced to: 

P,s{t) = X,{t)[l - XM, (16) 
since + Xs = 1. Furthermore, if the impact frequency is high, we get: 

At ^ dt dt ' ^ ^ 

Using eq. Uni and [T7| with eq. [T31we obtain the proportionality relation: 

at) oc ^ ^ (18) 

At this point we drop the lower indices e and 6 because they are no longer 
necessary. Equation describes the time evolution tribochemical activity 
to destroy (or create) the e phase. The same result applies to 6. 



3 Fe-Cu System 

We used the Fe55Cu45 system studied in a previous work P to perform an 
analysis of phase volume evolution as described in eq. [7| We used a Fritsch 
Pulverisete 5 planetary ball mill to perform two sets of experiments, with 

4 and 12 balls in the mill vial. The rotational speed used was 350 rpm for 
both experiments. The Fe-Cu samples were prepared by ball milling using 
commercial Fe and Cu metallic powders. In the course of milling, the iron 
atoms are incorporated into the fee copper matrix so that at any time during 
the process, there are fractions of the initial bcc Fe and fee Cu and an amount 
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of transformed fee FeCu solid solution. The observed trend is that of total 
incorporation of iron by the fee matrix resulting on fee FeCu alone jT| El • 
Therefore, in the crystallographic point of view, the system consists of a 
bcc phase (Fe) that decreases gradually and a fee phase (Cu + FeCu). We 
applied the model proposed in this contribution to analyze the bcc phase's 
kinetics of transformation. 

The grain size evolution presents a similar behavior for both bcc and 
fee phases and r^cc ~ '"/cc, (table [T] and fig. HJ. Thus, the evolution of 
tribochemical activity to destroy the bcc phase can be followed by (eq. ITH)) : 

ff>cc/,N dXhccjt) 1 , , 

^ ^' dt xum-xut)]' ^ ^ 

The evolution of phase concentration Xbcdt) can be deduced by the mag- 
netic hyperfine field distributions (BHF) obtained by Mossbauer spectroscopy 
(fig. 121) • Bcc concentration is the area under the spectra with BHF greater 
than 30 T. Figure El shows the evolution curves for the bcc phase volume in 
the course of milling. 

The tribochemical activities versus milling time for the two experiments 
described are shown in fig. jH The curves were calculated using eq. Uni and 
the data in Fig. El The scales, however, are not the same because impact 
frequencies, which are related to the number of balls, are different for each 
experiment. Figure El shows ^ vs. time for the two experiments rescaled by 
number of balls, in such a way that we get tribochemical activities per ball 
versus milling time per ball. 



4 Fe-N System 

The 2Fe + Fe2N — > Fe^N reaction obtained by MA studied in a previous 
work is also a good application for the model discussed in this contribution. 
The transformation was followed through the volume fraction of Fe2N since 
the existence domain of this paramagnetic phase is relatively narrow. In 
addition to this, the system 2Fe + Fe2N is two-phased, composed by a 
paramagnetic and a ferromagnetic phase; also the reaction product [Fe^N) 
is ferromagnetic. Therefore, the phase volume of the paramagnetic Fe2N can 
be easily separated by Mossbauer spectroscopy. It is the area of the doublet 
subspectra of Fig. El The Fe2N paramagnetic phase volume evolution is 
shown in Fig. [71 
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The Fe-N equilibrium phase diagram is composed by 4 crystallographic 
phases named a — Fe, 7' — Fe^N, e — Fe2.1-3.6N and 6 — Fe2N. At high 
temperature one also finds a gamma phase 7 — FexN and a epsilon nitrogen 
poor e — FciN. One can also find a metastable long range ordered phase a" — 
Fei6A^2- However, Fe-N phases obtained by milling are quite distinct. Milling 
of 7' — Fe^^N and 2Fe + Fe2N leads to a hexagonal phase which presents 
different nitrogen distribution. Moreover, milling of 4Fe + Fe^N leads to a 
hyperferrite instead of the well ordered a"-FeieN2 01111 11111201 1211 122 • The 
large domain of epsilon phases, the wide range of the Nitrogen distribution, 
and the nano crystalline character of these phases unable any attempt to 
distinguish nitrogen compositions in epsilon phases. What can be observed is 
that during most of the milling of 2Fe+Fe2N, the system is mainly composed 
of a — Fe, e — Fe2.1-3.QN and 6 — Fe2N . Previous studies about Curie 
temperature [211 aren't useful because differences in Nitrogen composition 
lead to diverse magnetic behavior |21]. Indeed, almost any variation on 
Nitrogen content leads to appearance of ferromagnetic phases, with only 
5 — Fe2N remaining paramagnetic. 

After just few minutes of alloying it is no longer possible to distinguish 
between different phases by X-ray diffraction. The local description provided 
by Mossbauer spectroscopy cannot distinguish Fe atoms without N as near 
neighbors from a — Fe. Even by transmission electron microscopy that sep- 
aration is no longer possible [21]. Therefore it is only possible to distinguish 
between ferromagnetic and paramagnetic iron in milled iron nitrides and on 
the magnetic point of view, the system is composed by two phases: one 
paramagnetic and one ferromagnetic. Since there are no direct link between 
tribochemical activity and crystallography (tribochemical activity is a pa- 
rameter to measure a given activity), this concept can be perfectly applied 
here. 

Either phase is suited for the kind of analysis performed in this paper. 
Due to experimental convenience, the paramagnetic phase was chosen. The 
evolution of tribochemical activity to destroy the paramagnetic phase can be 
followed by eq. UHl 

cP(+\ _ dXpit) 1 

^ " dt x,m-xp{t)]- 

The tribochemical activity versus milling time for the Fe-N system studied 
is showed in fig. |H1 The curve was calculated using eq. (201 and data in Fig. |3 
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5 Discussions 



The two curves in Fig. |3] are not in the same scale because the number of 
balls, which influences the number of impacts per unit time, is different in 
each case. We assume the number of impacts per unit time between balls 
and the vial wall equals the number of balls in the vial. In this way, there 
are at least n = 3 times as many impacts in the experiment with 12 balls. 
A correction to be made refers to the different number of impacts between 
balls in each experiment. It is virtually impossible to calculate such a number 
due to the nearly unpredictable dynamics of the balls in the vial. However, 
we empirically estimate this correction to be 1/5 of n. Therefore, there are 
about n = 3.6 as many impacts in the 12 balls case. To bring both curves 
to scale, we modify the 12 balls curve by multiplying time (and dividing ^) 
by n. Furthermore, we multiply the time axis (and divide the ^ axis) by 4 to 
obtain tribochemical activity vs. time, per ball. In this way, we get the two 
curves in Fig. El which represent the tribochemical activity time evolution we 
would obtain if there were only one ball in the vial. The curves in Fig. El are 
very similar indicating that the number of balls (impact frequency) does not 
play an important role on the final product. 

The behavior of the curves in Fig. El reveals four distinguished stages. 
The first one is a static activation stage where the system only stores energy 
and no chemical reaction occurs [T3|. The second stage is characterized 
by the start of reaction followed by an increase of reaction speed. In fact, 
the reaction occurs only in the bcc/fcc interface and so this acceleration is 
expected when total surface increases (grain size decreases) like in fig. ^ 
The third stage is characterized by a decrease in the reaction rate. During 
this stage, the amount of iron dissolved in the copper matrix approaches the 
solubility limit of ~60% which explains the decrease. The fourth 

and last stage is characterized by a zero ^ value revealing that stationary 
stage is reached. 

In contrast to what one can observe in the Fe-Cu formation process, the 
tribochemical activity for the Fe-N system shown in Fig. |H1 reveals large reac- 
tivity from the beginning of the process (^ is large when t = 0). Moreover, 
the reactivity seems to be proportional to the concentration gradient between 
the two phases {Fe and Fe2N) and the reaction rate increases due to an ever 
growing interface between Fe2N and Fe (grain size reduction). In fact, the 
Fe-Cu is an endothermic system and energy storage is necessary in order to 
trigger the reaction. This does not occur in the exothermic 2Fe + Fe2N re- 
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action. It can also be observed that the absolute value of ^pp is a monotonic 
time function showing that saturation has not been achieved, and that there 
was no reduction of reagents reactivity capacity. The reaction stops because 
of the reactants extinction. This is another contrast to the Fe-Cu system in 
which a stationary stage is reached. 

The presence of a stationary stage with zero ^ in the end of the Fe- 
Cu processing suggests us that there might be some unreacted amount of 
material left over. This fact calls for adjustment in eq. El Since the Fe-Cu 
reactions reached saturation, it is reasonable to believe that there should be 
a constant term in eq. El to account for the phase volume not transformed. 
Thus, the modified equation takes, the form: 

V{t) = ke-^'^ + R, (21) 

and plugging this expression into eq. [THl we obtain: 

Lndoit) « _^[2-fce''*^(l-i?)]" ^ ^ 

On the other hand, the Fe-N system never reached saturation indicating that 
there shouldn't be any unreacted material left over and therefore R should 
be zero for exothermic reactions reducing eq. |^ to the original one (eq. E)) 
and eq. I^to: 

C.oit) oc -YZTk^- (23) 

It is useful to analyze the limiting behavior of eqs. EH and ESI For large 
t, eq. E21 becomes: 

~ canst. + {R-R^)e^^- ^^^^ 

For endothermic processes, R > and |^| ~ 1/e*'' — * 0, as expected. In 
eq. ESI the numerator dominates and |^| ~ t^^^ oo, which is consistent 
with exothermic reactions. For small values of t, the denominator in eq. E21 
and ESI is constant and then the numerator is dominant. The time derivative 
of ^ when t is small goes like — Therefore, for c > 2, dC,/dt is zero in the 
beginning of the reaction, which is congruent with an energy storage stage 
which should be present in endothermic processes. On the other hand, for 
c < 2, \dC,/dt\ is maximum when t —>■ which is compatible with exothermic 
reactions. 
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To test these limits, fittings of tlie tribochemical activity time evolution 
curves in Figs. and |H1 to eqs. |221 and 1221 respectively, were performed. It 
is important to point out that eqs. |221 and 1221 are result of an empirical 
assumption and there is no argument to sustain that the fittings made here 
are, by any means, the best ones. One could obtain fittings as good as 
these ones using many different expressions. However, although the aim 
of the analysis made here is not to establish the absolute validity of the 
law in eqs. |21 and |^ it provides further insights on the usefulness of the 
tribochemical activity approach. 

Table |21 shows the values of coefficients k, b, c, and R obtained from 
numerical fitting of tribochemical activity time evolution curves to eqs. |221 
and 1221 where the coefficient R applies to the former but not to the latter. 
Comparing the values obtained from the two Fe-Cu experiments (with 4 and 
12 balls) we see that the values of coefficient b are quite different, which is 
in accordance to the proposed relation of them to "work conditions" . The 
values of c from both experiments are very similar. This was expected since c 
should be connected to impact energies and both experiments were performed 
with rotational frequency of 350 rpm. As a consequence, approximately the 
same amount of energy is transferred to the powder in both experiments and 
the local kinetics of both reactions are the same, despite the different number 
of balls. This analysis strengthens the suggestion that the shock frequency 
is only a milling condition and does not play a role in the energy transfer 
process, i.e., it is not a fundamental parameter of the process. The values of 
k and R are virtually the same since one process differs from the other by 
only the time scale. Moreover, according to tab. |21 c is less than 2 for the 
exothermic Fe-N system in contrast to the endothermic Fe-Cu system where c 
is greater than 2. These numbers are in accordance to the previous discussion. 
In this way, not only is c connected to impact energies as suggested in [T2|, 
but also reflects the endothermic or exothermic nature of the reaction. 

6 Conclusions 

The formulation of the tribochemical activity concept provides a clear way 
to understand the evolution of phase transformations during MA. This is an 
initial attempt to describe the growth laws from global kinetics. Using this 
description it is possible to follow grain sizes, gradient of concentration and 
milling condition without chemical conditions. 
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The time evolution of tribochemical activity allows the identification of 
four distinct stages in the milling of Fe-Cu system: energy storage, reaction 
initiation and rate increase, reaction rate decrease, and stationary. These 
stages are in perfect agreement to the experimental observations. The energy 
storage stage is not identified in the 2Fe + Fe2N — > Fe^N reaction because 
of the exothermic nature of this system. Moreover, the monotonous character 
of the ^{t) function shows that saturation is not achieved. 

The striking difference between ^ behavior for Fe-N and Fe-Cu reflects 
the respective exothermic and endothermic nature of the reactions. Achieve- 
ment of a saturation stage in Fe-Cu processes suggests the existence of some 
amount of unreacted material left over. This fact calls for modification on 
the global kinetics law proposed in [T2]. The small adjustment introduced 
extends the validity of the law to endothermic and exothermic processes and 
imposes a threshold to the characteristic exponent. 

Further analysis of the two Fe-Cu tribochemical activity curves suggests 
that impact frequency does not play an important role in the final product 
process. This is stressed by the values obtained from the numerical fitting of 
the phase volume and ^ time evolutions. 

At this point, it is really hard, and certainly premature, to assert that the 
tribochemical activity approach is definitive, but it definitely makes it clearer 
that the right path to fully understand the MA processes should bend more 
and more toward its description in terms of growth laws. Such description 
provides some clues about which parameters actually drive MA mechanisms, 
and also provides a means to separate out surface, mechanical impacts, and 
thermodynamics terms. The rather different tribochemical activity behav- 
ior for the Fe-Cu and Fe-N systems might explain the broad diversity and 
incompatibility among various MA models present in the literature. 
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Table 1: Grain size for some Fe-Cu milled samples with 12 balls and 350 
rpm. This results were obtained refining the x-rays patterns using Rietveld 
method. 



time(h) 


grain size (nm) 


bcc phase 


fee phase 


1 


74 


44 


2 


41 


38 


4 


33 


24 


8 


15 


14 


12 


14 


14 


16 




13 



Table 2: Coefficients obtained by fitting Figs. |31 and |H1 to eqs. |221 and ESI 
respectively. 





k 


b 


e 


R 


Fe-Cu (4 b) 


0.479 


0.003 


2.729 


0.075 


Fe-Cu (12 b) 


0.485 


0.121 


2.658 


0.069 


Fe-N 


0.497 


0.106 


1.661 
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Figure 1: Grain size evolution of bcc and fee phases during milling for Fe- 
Cu milled with 350 rpm and 4 balls revealing a similar behavior for the two 
phases. 
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(a) 350 rpm - 4 balls 
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(b) 350 rpm - 12 balls 



Figure 2: BHF distributions for some Fe-Cu milled samples showing a clear 
separation between Fe-fcc and Fe-bcc phases (low and high BHF values re- 
spectively) 
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Figure 3: Bcc relative volume phase for the Fe-Cu system obtained by 
Mossbauer spectroscopy. 
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Figure 4: Tribochemical activity to destroy the bcc phase in the Fe-Cu system 
as a function of milhng time. The curves were calculated from data in Fig. El 
using eq. [THl 
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Figure 5: Tribochemical activity to destroy the bcc phase in the Fe-Cu system 
as a function of miUing time rescaled by the number of balls. 
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Figure 6: Mossbauer spectra at room temperature for milled Fe2-/V + 2Fe 
up to (a) 5min; (b) 20min; (c) Ih; (d) 2h; (e) 4h and (f ) 16h 
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Figure 7: Paramagnetic relative volume phase for the Fe-N system obtained 
by Mossbauer spectroscopy 
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Figure 8: Time evolution of tribochemical activity to destroy the paramag- 
netic phase in the Fe-N system. The curve was calculated from data in Fig. [7| 
using eq. 1201 
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